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Using backwards ray tracing, we study the shadows of Kerr black holes with scalar hair (KBHsSH). 
KBHsSH interpolate continuously between Kerr BHs and boson stars (BSs), so we start by inves¬ 
tigating the lensing of light due to BSs. Moving from the weak to the strong gravity region, BSs 
- which by themselves have no shadows - are classified, according to the lensing produced, as: (i) 
non-compact, which yield no multiple images; (ii) compact, which produce an increasing number of 
Einstein rings and multiple images of the whole celestial sphere; (iii) ultra-compact, which possess 
light rings, yielding an infinite number of images with (we conjecture) a self-similar structure. The 
shadows of KBHsSH, for Kerr-like horizons and non-compact BS-like hair, are analogous to, but 
distinguishable from, those of comparable Kerr BHs. But for non-Kerr-like horizons and ultra¬ 
compact BS-like hair, the shadows of KBHsSH are drastically different: novel shapes arise, sizes 
are considerably smaller and multiple shadows of a single BH become possible. Thus, KBHsSH 
provide quantitatively and qualitatively new templates for ongoing (and future) very large baseline 
interferometry (VLBI) observations of BH shadows, such as those of the Event Horizon Telescope. 

BAGS numbers: 04.20.-q, 04.70.Bw 04.80.Cc 


Introduction, 100 years after General Relativity (GR) 
was formulated, we finally face a realistic prospect of test¬ 
ing one of its most dramatic consequences: black holes 
(BHs). The evidence for astrophysical BHs, gathered for 
over half a century, has built a strong case [T] , but it could 
not confirm the existence of event horizons, the defining 
property of BHs. The near future promises to open up 
a new channel of observation - gravitational waves - and 
deliver electromagnetic measurements of unprecedented 
precision, hopefully clarifying this central issue [2] . 

A particularly exciting prospect is the use of VLBI 
techniques to resolve the angular scale of the event hori¬ 
zon for some supermassive BH candidates and determine 
the corresponding BH shadow 0 SI. Its observation 
would probe the spacetime geometry in the vicinity of the 
horizon and consequently test the existence and proper¬ 
ties of the latter [5]. It is therefore timely to study BH 
models that yield phenomenological deviations from the 
paradigmatic GR BH, described by the Kerr metric. 

One approach is to parameterize families of metric de¬ 
viations from Kerr [Mil]. Another approach is to use 
exact solutions of GR (or generalizations thereof) yield¬ 
ing deviations from Kerr (e.g. HlHIll). Exact solutions 
with physically reasonable and astrophysically plausible 
matter sources, however, are scarce; but Kerr BHs with 
scalar hair (KBHsSH) [15] are arguably one such model. 
These are exact solutions of Einstein’s gravity minimally 
coupled to a massive complex scalar field, and interpolate 
between Kerr BHs and gravitating solitons - boson stars 
(BSs) US] - suggested as dark matter candidates (in the 
Newtonian limit) and BH mimickers [TTHTO] . 

In this letter we show that the shadows of KBHsSH are 
distinguishable, or even drastically different, from those 
of Kerr BHs, and can thus yield new templates for the 
ongoing and future VLBI searches [20] of BH shadows. 


The solutions, KBHsSH can be expressed by a station¬ 
ary and axi-symmetric line element, in spheroidal coordi¬ 
nates (t, r, 0, (j)) [mUl], together with the (mass /i) scalar 
field At = (j){r, 0) ^ where w is the frequency and 

m G z+ is the azimuthal harmonic index. The metric 
functions and (p are determined numerically by solving 
five coupled, non-linear PDFs [21]. For m ^ 1, both BSs 
and KBHsSH solutions can be obtained. For m = 0, 
spherical BSs exist. The space of solutions for m = 0,1 
is summarized in Fig. BSs exist for a limited range of 
w along spiraling curves. KBHsSH exist inside an open 
set, bounded by the m = 1 BS curve, a set of Kerr BHs 
and the set of extremal KBHsSH. 
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FIG. 1: BS solutions with m = 0,1 (red dashed and solid 
lines) in an ADM mass Madm vs . scalar field frequency w 
diagram (in units of /a). For m = 1, KBHsSH exist within 
the shaded region. Points 1-12 (I-HI) correspond to the BSs 
(KBHsSH) exhibited below. LR=liglit ring. 
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Setup, Our goal is to compute photon geodesics in 
the above geometries and obtain, at a given observation 
point, the distorted apparent sky, due to the gravitational 
leasing of BSs, as well as the shadow, when a horizon is 
present. To interpret the patterns obtained, we divide 
the “celestial sphere” light source into four quadrants, 
each painted with a different color. On top of these, a 
grid of constant longitude and latitude (black) lines is 
introduced, with adjacent lines separated by 10° - Fig.|^ 
(left). This setup mimicks closely the one in [22] . 




FIG. 2: (Left) The full celestial sphere. (Right) The viewing 
angle of the off-centred observer {O). 


The observer, henceforth denoted as O, is placed in¬ 
side the celestial sphere at an off-centered position to be 
specified below, and its viewing angle captures ^ 100° 
of that sphere - Fig. (right). The point on the celes¬ 
tial sphere immediately in front of O (dubbed F) lies at 
the intersection of the four colored quadrants and it is 
marked with a white dot. 

From the O’s position we span its viewing angle 
by performing, numerically, backwards ray tracing of 
2000 X 2000 photon trajectories. We integrate these null 
geodesics until they either reach a point on the celestial 
sphere or they hit the BH horizon (when it is present). 
The latter directions form the shadow [29]. O’s image 
is obtained upon a projection onto its local frame, by a 
method similar to that in [22] . 

Quantitative shadow parameters. To analyse the 
shadows obtained below we introduce six parameters, 
{Dc , Dx , Dy cTKerr}, mostly following [23| . Let 

O’s image be parametrized by the Cartesian coordinates 
(x,^), obtained from multiplying the observation angles 
by the circumferential radius R [30| at O’s position. The 
origin of this coordinate system, O, points at the (un- 
lensed) white dot on the celestial sphere. The centre 
of the shadow, C, has xc = + ^min)/2, where 

Xmax and Xmin (^max and ^min) are respectively the maxi¬ 
mum and minimum abscissae (ordinates) of the shadow’s 
edge. Due to reflection symmetry for observations on the 
equatorial plane, = 0. C and O need not coincide; 
the displacement^ Dc = \xc\^ measures their difference. 
The width and height of the shadow are, respectively, 
Dx = ^max ^min and Dy — ^max ^min • 

A generic point P on the shadow’s edge is at a distance 


r = {yp‘^ + {xp — xc)‘^y^‘^ from C. Let a be the angle 
between the line CP and the x axis; the average radius 


is r = r(a) da/27r and the deviation from sphericity 

nl/2 


IS ar = 


C (^(“) “ 'f 


. Finally, 


(JK = 



r{a) - rKerr(<^) 

“^Kerr (<^) 


2 

da , 


( 1 ) 


is the relative deviation from a comparable Kerr BH, 
either with the same ADM mass and angular momen¬ 
tum, Madm, ot with the same horizon quantities, 

Mh, Jr (as long as the Kerr bound is not violated for the 
comparable Kerr BH). These comparable Kerr BHs are 
denoted below as Kerr=ADM and Kerr=H, respectively. 
C is made to coincide for both BHs in Q. 

Leasing by spherical BSs. We first look at the leasing 
due to spherically symmetric BSs - Fig. We set O on 
the equator and always at = 22.hj\i, for the different 
BS solutions. Throughout we take G = 1 = c and in the 
following, unless otherwise stated, y = 1. 

Starting from vacuum {i.e. w = 1) we find a set 
of non-compact BSs, for which any meaningful effec¬ 
tive radius is large as compared to the corresponding 
Schwarzschild radius. These are “weak gravity” solu¬ 
tions and the corresponding leasing is illustrated by the 
BS with = 0.95 [31] in Fig. where only a small 
distortion of the background is observed. Moving further 
along the spiral, an Einstein ring appears at ^ 0.94. 

The Einstein ring is formed by the leasing of F, and it 
encloses two inverted copies of a region around F. This is 
illustrated in Eig. [fusing a BS with = 0.9. The ap¬ 
pearance of the first Einstein ring defines the transition 
from non-compact to compact BSs. 



FIG. 3: Leasing by spherical BSs. From left to right: (top) 
'^ 1 ^ 2,3 = 0.95; 0.9; 0.8; (bottom) — 0.8; 0.84 (zoomed). 
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Moving further along the spiral the region duplicated 
inside the Einstein ring becomes an increasingly larger 
part of the celestial sphere, as can be observed for the 
BS with = 0.8 in Fig. ^ Shortly after the spi¬ 
ral’s first backbending, at ^ 0.767 (c/. Fig. 

the full celestial sphere becomes duplicated, starting at 
the BS with — 0-77. At, and beyond, this point, 

two further Einstein rings emerge - now corresponding 
to lensings of the point immediately behind the observer 
(dubbed B). This is illustrated by the bottom left panel 
in Fig. for a BS with = 0.8. In between these 

two new Einstein rings, further pairs of Einstein rings can 
be seen to appear, progressively, further inside the spiral. 
The first such pair is illustrated in Fig. [^for = 0.84. 
Each new pair of Einstein rings corresponds to images of 
either F or 5, in an alternating fashion, and to a further 
complete copy of the full celestial sphere. 

An infinite number of copies, and a corresponding self¬ 
similar structure, is expected to arise when a light-ring 
- corresponding in this case to a photon sphere - ap¬ 
pears [32] , marking the transition from compact to ultra¬ 
compact m BSs. This occurs well inside the spiral, on 
the third branch (after the second backbending), start¬ 
ing at the BS solution with = 0.842, marked as the 
blue star ‘LF’ point on the m = 0 spiral of Fig. 

Lensing by rotating BSs, We now turn to rotating 
BSs (spin axis pointing up). Starting from vacuum, we 
again find a region of non-compact BSs, i.e without mul¬ 
tiple images. Two differences, however, with respect to 
the top left panel of Fig. are: an asymmetric lensing, 
with an amplification of the side rotating away from O 
and the slight shift of point F to the left, due to frame 
dragging; these are illustrated in Fig. for a BS with 
= 0.95. At 0.92 an Einstein ring appears, 

startingthe set of compact rotating BSs. This is the case, 
in Fig. H for the BSs with = 0.85; 0.8. The ring en¬ 
closes again two inverted copies of part of the celestial 
sphere, but it is now elliptic and the duplicated image 
of the side rotating towards the observer is suppressed. 
The inversion shifts point F to the right. 

Further following the spiral, new Einstein ‘rings’ ap¬ 
pear, just as in the case of spherical BSs, except that 
instead of an ‘0-shape’, they have a squashed ‘D-shape’. 
An example, for a BS with = 0.75 is shown in Fig.|^ 
Then, a light ring on the equatorial plane appears at 
0.747. Beyond this point multiple (presumably 
infinitely many) images of the celestial sphere arise, with, 
we conjecture, a fractal structure. This is illustrated W 
the BSs with = 0.7; 0.65 and = 0.7 in Fig.j^ 

Shadows of KBHsSH. A KBHSH may be regarded 
as a BS around (and co-rotating with) a central hori¬ 
zon. The latter may be non-Kerr-like^ i.e violate the 
Kerr bound in terms of horizon quantities [25] and the 
former may have strong lensing effects, as seen above. 
Consequently, it is expectable that KBHsSH with ultra¬ 
compact BS-like hair and with a non-Kerr-like central 



FIG. 4: Lensing by rotating BSs with m — 1. From left to 
right: (top) g = 0.95; 0.85; 0.8; (middle) = 0.75; 0.7; 
(bottom): = 0.65; = 0.7. 


horizon will have unfamiliar shadows. This expectation 
is confirmed in Fig. [^ where we exhibit the shadows and 
lensing for three examples of KBHsSH (together with 
four ‘transition’ examples) and for comparable Kerr BHs. 
The three main examples are configuration I-HI in Fig.[^ 
and their physical quantities, in units of g when dimen¬ 
sionful, are summarized as follows: 



Madm 

Mh 

-^ADM 

Jr 

Mh 


-^ADM 

Jh 


Madm 

-^ADM 

^ADM 


I 

0.415 

0.393 

0.172 

0.150 

95% 

87% 

0.999 

0.971 

H 

0.933 

0.234 

0.740 

0.115 

25% 

15% 

0.850 

2.10 

HI 

0.975 

0.018 

0.85 

0.002 

1.8% 

2.4% 

0.894 

6.20 


In Fig.[^ we have scaled g in each case so that Madm 
is the same for all KBHsSH presented, and O sits at 
R = 15Madm- Comparing the shadows in this way bears 
more significance for real observations, where Madm of 
the BH is fixed by the data. 

The top left panel of Fig. shows the shadow of con¬ 
figuration I (shadow I for short). For this KBHSH the 
‘hair’ contains only 5% of the total mass; the ‘nearby’ 
BSs in Fig. are non-compact; and the central horizon 
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FIG. 5: From left to right: shadows of (top) configuration 
I, Kerr=ADM and Kerr=H; row) conhguration II and 

Kerr=ADM; (3^^^ row) transition conhgurations between II and 
III (detail); (bottom) conhguration III and Kerr=ADM. 


is Kerr-like (Jh/M^ < 1). The shadow, albeit qualita¬ 
tively familiar, is nevertheless distinguishable from that 
of the Kerr=ADM BH, which is exhibited in the top mid¬ 
dle panel. The latter is slightly larger and more T)-like - a 
characteristic of extremal Kerr BHs. Shadow I turns out 
to be closer to the one of the Kerr=H BH, exhibited in the 
top right panel. This observation can be quantitatively 
checked: ctk = 4.81%/0.52%, taking Kerr=ADM/Kerr=H 
for the comparable BH, cf. Table 1. 

New types of BH shadows, quite distinct from those of 
Kerr=ADM BHs, appear on the left 2^^ and 4^^ row pan¬ 
els of Fig. corresponding to shadow H and HI. In both 
cases, the central BH is non-Kerr-like (Jbh/^bh ^ ^)- 
Shadow H is smaller (average radius of ^75%) than that 
of the Kerr=ADM BH. It is also more ‘square’, with a 
larger normalized deviation from sphericity. Shadow HI 
is remarkably distinct. Its central BH has Ju/M^ ^ 6, 




D (j Dx d^y ^ 

^(%) 

<Jk{%) 

Shadow I 


2.07 8.48 9.33 4.48 0.170 

3.8 

4.81/0.52 

Kerr=ADM 


2.30 8.73 9.76 4.67 0.217 

4.64 

1.21 

Kerr=ADM 

A 2.38 8.66 9.86 4.70 0.260 

5.54 

0 

Kerr=H A 


2.07 8.48 9.36 4.50 0.180 

3.99 

0 

Shadow 11 


2.39 7.14 6.93 3.60 0.118 

3.29 

25.5 

Kerr=ADM 

A 

1.79 9.32 9.86 4.82 0.103 

2.15 

0 

Shadow HI 

1.79 5.30 4.67 1.63 0.838 

51.3 

68.1 

Kerr=ADM 

A 

1.92 9.22 9.86 4.80 0.125 

2.60 

0 


TABLE I: Parameters for Kerr shadows with ‘A’ are com¬ 
puted from the analytic solution [3]; ctr is always computed 
with respect to such solution. The second line in the table, 
computed for a Kerr BH generated numerically and using the 
same ray tracing code as for KBHsSH, estimates the numeri¬ 
cal error (~ 1%) of the KBHsSH shadows. 


allowed by the ‘heavy’ hair that it is dragging {cf. the dis¬ 
cussion in [25] ). The leasing of this hair resembles closely 
that of the ultra-compact BS on the bottom left panel of 
Fig. a Interestingly, multiple disconnected shadows of 
the (single) BH occur: the largest ones (besides the main 
‘hammer-like’ shadow) are two eyebrows [26l|27|, at sym¬ 
metric positions above and below the main shadow; but 
we have detected many other smaller shadows, hintiM 
again at a self-similar structure. On the 3^*^ row of Fig.[^ 
the shadows of four solutions in between configurations H 
and HI illustrate the transition between them. Finally, 
we remark that the shadows of KBHsSH can have ar¬ 
bitrarily small sizes by considering solutions arbitrarily 
close to the BS curve in Fig. 

Remarks, KBHsSH can lead to qualitatively novel types 
of shadows in GR, as shown by shadows H and HI. Even 
for KBHsSH close to Kerr, their shadows are distinguish¬ 
able from the latter, with the same asymptotic quantities, 
as illustrated by shadow I. Regardless of the astrophysi- 
cal relevance of these solutions - which is unclear - they 
can yield new templates with small or large deviations 
from the Kerr shadows, hopefully of use for VLBI ob¬ 
servations. An exhaustive analysis of KBHsSH shadows 
spanning the space of solutions in Fig. and at differ¬ 
ent observation angles, will be presented elsewhere, for 
producing such templates [28|. But the examples herein 
already raise a challenge to the parameterizations of de¬ 
viations from Kerr suggested in the literature I2HII] : can 
they describe shadows with such large deviations? 

Besides the peculiar shape of some of the shadows ex¬ 
hibited, this model has one general prediction: smaller 
observed shadows than those expected for Kerr BHs with 
the same asymptotic charges. Indeed, a ‘smaller’ central 
BH seems a natural consequence of the existence of hair, 
carrying part of the total energy. 

Finally, for the setup herein, the redshift, which de¬ 
pends only on the source’s and O’s positions, is constant 
throughout the image and has been neglected. 
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